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We perform a sensitivity study of dark energy constraints from galaxy cluster surveys to uncer- 
tainties in the halo mass function, bias and the mass-observable relation. For a set of idealized 
surveys, we evaluate cosmological constraints as priors on sixteen nuisance parameters in the halo 
modeling are varied. We find that surveys with a higher mass limit are more sensitive to mass- 
observable uncertainties while surveys with low mass limits that probe more of the mass function 
shape and evolution are more sensitive to mass function errors. We examine the correlations among 
nuisance and cosmological parameters. Mass function parameters are strongly positively (nega- 
tively) correlated with ^de {w). For the mass-observable parameters, Qde is most sensitive to 
the normalization and its redshift evolution while w is more sensitive to redshift evolution in the 
variance. While survey performance is limited mainly by mass-observable uncertainties, the current 
level of mass function error is responsible for up to a factor of two degradation in ideal cosmological 
constraints. For surveys that probe to low masses {10^^'^h~^ Mq), even percent-level constraints 
on model nuisance parameters result in a degradation of ~ \/2 (2) on f^DE (w) relative to perfect 
knowledge. 



I. INTRODUCTION 

The spatial abundance of galaxy clusters is a poten- 
tially powerful approach to test the nature of dark en- 
ergy [IHZ]- As is the case for all dark energy probes, 
modeling and other systematic uncertainties in the data 
analysis must be identified and controlled. The primary 
uncertainties for cluster counts involve theoretical mod- 
eling of the halo space density as a function of mass — 
the mass function — along with uncertainties associated 
with modeling the halo mass selection function of a spe- 
cific survey. 

While recent studies have derived consistent and com- 
petitive cosmological constraints using X-ray [SHIO] and 
optical [m [12] cluster surveys, the larger volume and 
improved sensitivity of upcoming surveys motivate a 
stronger effort to study sources of systematic error. Early 
work on forecasting cluster constraints focused attention 
on systematics arising from the mass-observable relation 
([T3HI1]: etc.). These errors fall into two categories: i) 
bias arising from uncertain knowledge of the normaliza- 
tion and ii) mass variance at fixed observable signal or, 
equivalently, scatter in the mass-signal relation. The red- 
shift evolution of these effects are particularly important 
[I6l - fl9| . If the scatter is large, the functional form of the 
mass-observable likelihood becomes important ^0] . 

Uncertainties in the mass-observable relation can be 
characterized with a combination of empirical and the- 
oretical approaches. For example, the slope and nor- 
malization of the intracluster gas thermal energy can be 
probed with X-ray and Sunyaev-Zel'dovich (SZ) observa- 
tions pTH27] while the variance and evolution with red- 
shift can be predicted in a model-dependent manner us- 
ing multi-fluid simulations p5H55] . For example, a 'pre- 



heated' gas treatment that matches the low-redshift X- 
ray luminosity-temperature relation predicts weakly in- 
creasing scatter in the integrated thermal SZ signal at 
fixed mass, from 13% at 2; = to 18% at z = 1 [55] . 

Theoretical uncertainty in the mass function has not 
been extensively studied in previous Fisher forecasts. 
This is partly because mass-observable uncertainties are 
typically dominant in comparison, and partly because 
mass function calibrations by the simulation community 
have been evolving. 

The original mass function of Press & Schechter |37j . 
derived with support of 1000-particle N-body simula- 
tions, has evolved into a number of forms with billion- 
particle simulation support, among them the Sheth- 
Tormen |5S], Jenkins [39;, Evrard [40 , Warren gT] and 
Tinker [42] mass functions. The Warren parameteriza- 
tion has received a recent update based on an ensemble 
of large- volume simulations^. Correction for systematics 
associated with initial condition generation and other ef- 
fects lead to an upward revision in the mass function of 
> 10 percent at mass scales probed by SZ surveys 

Different measures of mass, based on either particle 
percolation or spherical filtering, are employed by the 
aforementioned studies (see piHi^ for discussion of these 
mass measures) , but all use the filtered linear power spec- 
trum, a(M), as a similarity variable for expressing halo 
counts and clustering. While the Tinker calibration [15] 
improved the statistical accuracy of the mass function to 
the 5% level in number density, the effects of gas dynam- 
ics, absent in that work, may produce deviations that 
are larger than this [ITj- We consider ~ 10% in space 
density as a reasonable estimate of the current level of 
uncertainty in the mass function at cluster scales. 

The cited mass function calibrations assume a stan- 
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dard ACDM cosmology, with a Gaussian initial density 
field evolved under general relativity with a cosmological 
constant. Extensions to non-standard assumptions are 
emerging (see e.g. [JBHSU]), but we do not consider them 
here. 

In this paper, we extend previous Fisher studies of dark 
energy constraints by including theoretical uncertainties 
in the halo mass function and clustering bias and by em- 
ploying a more generous treatment of bias and scatter 
evolution in the mass-observable relation. We frame the 
analysis in terms of two idealized surveys roughly pat- 
terned after the South Pole Telescope^ (SPT) SZ survey 
and the Dark Energy Survey'^ (DES) optical survey. Our 
baseline idealizations, which assume perfect halo selec- 
tion (100% completeness) above a redshift- independent 
mass threshold and constant mass-observable scatter, 
provide a simple, rational basis upon which complexi- 
ties, in the form of nuisance parameters discussed below, 
can be added. The "SZ-motivated" survey has a high 
mass threshold and small mass-observable scatter while 
the "optical-motivated" case has a lower mass threshold 
and larger mass scatter. We refer to these simply as 'SZ' 
and 'optical' for the remainder of the paper. 

Our main aim is to study performance degradation in 
f2DE and w constraints as a function of prior parameter 
knowledge. We also provide a first look at parameter 
correlations for the full model, finding them (unsurpris- 
ingly) complex. While writing up this work, we learned 
of a similar study by Wu, Zentner & Wechsler [51], who 
address the effect of mass function systematics on the 
dark energy figure of merit from cluster surveys. Their 
approach is complementary to ours, in that they adopt a 
set of independent nuisance parameters to model system- 
atic error in binned representations of the mass function 
and halo bias. While this non-parametric approach al- 
lows them to more easily identify regions of mass and 
redshift space that most sensitively affect dark energy 
constraints, a disadvantage is that it disregards correla- 
tions between different mass/redshift bins. Wu et al. also 
explore a time-dependent equation of state, ^(a), and 
express their results in terms of a figure of merit degra- 
dation. They do not emphasize the interplay between 
mass-observable and mass function systematic errors, as 
we do here. 

The paper is organized as follows. In Sec. |lT]we briefly 
review the formalism for extracting dark energy con- 
straints from cluster counts and variance in counts and 
present our parameterization of the mass-observable re- 
lations, the mass function and galaxy bias. We present 



results in Sec. Ill offer a critique in Sec. IV and conclude 
in Sec. |Vl 



^ http://pole.uchicago.edu/ 
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II. DE FROM CLUSTER COUNTS AND 
CLUSTERING 

The subject of deriving cosmological constraints from 
cluster number counts and clustering of clusters has been 
treated extensively in the literature (see e.g. [TSHIS]). In 
this section we follow closely the approach described in 

m- 

The comoving number density of clusters at a given 
redshift z with observable in the range M"^^^ < Mobs < 
^o"bt' is given by 

/-^^o-'dMobs fdM dn 



Mobs J M dluM' 



Uncertainties in the redshifts distort the volume ele- 
ment. Assuming photometric techniques are used to de- 
termine the redshifts of the clusters (hereafter photo-z's) 
and a perfect angular selection the mean number of clus- 
ters in a photo-z bin < < zf^-^ is 



'*'dzP / dz^n„H/*^(17)p(zP|z) (2) 
./ dz 



where Wf^{Vl) is an angular top hat window function. We 
parameterize the probability of measuring a photometric 
redshift, zP, given the true cluster redshift z as [18] 



p(zP|z) 



1 



exp [-yHz^]] 



(3) 



where 



-.bias 



V2a, 



(4) 



^.bias jg ^YiQ photometric redshift bias and az is the scatter 
in the photo-z's. We fix the photo-z bias and scatter at 
and 0.02 throughout this paper. 

The sample covariance of counts m^.i is, given by |52| 



d^k 



(5) 



(27r) 



.W*{k)W,{k)^P,{k)P,{k), (6) 



where ha^i{z) is the average cluster linear bias defined as 

&Q,i(z) 



1 /"^^"bt' dMobs f dM 



na,i{z) Jm^^_ Mobs J M 

dna,i{z) 



dlnM 



-6(M;z)p(Afobs|M). (7) 



W*(k) is the Fourier transform of the top-hat window 
function and Pi{k) is the linear power spectrum at the 
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centroid of redshift 
b{M;z) 
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bin i. We present our choice for 
when we discuss the parameterization 
of the errors in the mass function and galaxy bias. We 
only calculate covariance terms for which i ^ j since off- 
diagonal terms are negligible. 

Following |18j . we find that the window function 
W* (k) in the presence of photo-z errors is given by 



2 exp 



ikn 



H, 



exp 



2Hf 



sin(fc||fr,/2) Ji{kx_ri6s) 



(8) 



Here 



r{zf ) is the angular diameter distance to the 



photo-z bin, and Sr. 



r{zf). Similarly, 



(T^(z). We assumed that H{z), z^'''''{z), 



H{z), 



= i7(zf) 

and az{z) are constant inside each bin. The variables fc|| 
and fc_L represent parallel and perpendicular components 
of the wavenumber k relative to the line of sight. 
Define the covariance matrix of halo counts 



Cij — Sij + TfiiSij 



(9) 

where fhi is the vector of mean counts defined in Eq. ([2| 
and Sij is the sample covariance defined in Eq. ([6]). The 
indices i and j here run over all mass and redshift bins. 
Assuming Poisson noise and sample variance are the only 
sources of noise, the Fisher matrix is, [THl [Ml [M] 



i,^-H-Tr[C-iS,„C-iS,^], (10) 



where the "," denote derivatives with respect to the 
model parameters. The first term on the right-hand side 
contains the "information" from the mean counts, to. 
The Sij matrix only contributes noise to this term, and 
hence only reduces its information content. The second 
term contains the information from the sample covari- 
ance. 



A. Systematics in the mass— observable relation 

We introduce six degrees of freedom in the sector of the 
model that links halo mass to an observable signal. We 
assume a log-normal form for the probability of measur- 
ing an observable signal, denoted Mpbs, given true mass 
M, 



p(Mobs|M) = 



M 



exp [-a;2 (Mobs)] , (11) 



where 



, lnMobs-lnAf-lnMbias(Mobs,z) 
x(Mobs) = -7= — ^ . (12) 



We model systematic error in the mass proxy by intro- 
ducing a redshift-dependent bias and variance 



lnA/bias(^) 



Bo + Si(l + z), 



^In Af 



(z) = al+Y,Srz\ 



(13) 
(14) 



i=l 



where Bq, Bi, aa and the three variance coefficients are 
assumed to be independent of mass. Our default assump- 
tion is that ctq is non-zero, with values discussed in Sec. 



II D while the remaining parameters have a fiducial value 
of zero. All six parameters are taken as degrees of free- 
dom and varied in the Fisher analysis presented in Sec. 

ttni 



B. Systematics in the halo mass function and bias 

In the sector of the model describing the halo mass 
function, we add ten more degrees of freedom. We write 
the space density of halos as 



dn 
dM 



rflncr 



M dM 

and adopt the Tinker parameterization of /(cr) |42j 



/(^) - A 





— a 




[(9 




hi 









(15) 



(16) 



Following [42], we allow the first three parameters of f{a) 
to vary with redshift, so that 



^(z) = Ao(l + z)^^ 
a{z) = ao(l + 2)°=" 
6(z) -5o(l + z)-° 



(17) 
(18) 
(19) 



For fiducial parameters, we adopt the values of [42] 
at A = 200: Aq = 0.186, = -0.14, ao = 1.47, 
a, = -0.06, bo = 2.57, logio(a) = ( ;JaA5) )'-'> and 
c = 1.19. As Tinker et al. [35] explain, A controls the 
overall amplitude of /(cr), a controls the tilt, and b sets 
the mass scale where the power law in /(cr) becomes sig- 
nificant. 

We adopt the b{M, z) fit of |38| for the galaxy bias 



6(M, z) 



1 



1 



2pc 



(20) 



and choose the fiducial values for the parameters to be 
ac = 0.75, 5c = 1.69, and pc = 0.3. 

In total, the mass function and bias introduce ten addi- 
tional parameters, and we consider all of these as degrees 
of freedom in the Fisher analysis presented in § |III| The 
assumption that the bias is independent of the mass func- 
tion is very conservative. Manera et al. [55] show that, 
in the range of scales we are interested in, the bias can 
be predicted to roughly ^ 10% accuracy given the mass 
function. 
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TABLE I: Halo modeling nuisance parameters 



Class 


Name Fid. Value 


Notes 






u.u 


constant bias 




til 


u.u 


z-dependent bias 


Mobs 


„2 

f^o 


U.U4,U.DZ0 


constant variance 




c, 
oi 


0.0 


variance 




02 


0.0 


z'^ variance 




<J3 


0.0 


z"^ variance 




^0 


0.186 


z- independent 




ao 


0.147 


;5 


Mass 




2.57 


)J 


Function 


c 


1.19 








-0.14 A{z) = ^0(1 + 2)-*=" 




Ox 


-0.06 


a{z) = ao(l + 2)°- 




a 


0.0107 


= bo(l + 2)" 






0.75 


z- independent 


Bias 


5c 


1.69 








0.30 


55 



C. Fiducial Parameter Values 

The fiducial values of the sixteen nuisance parameters 
are summarized in Table |T] These parameters control the 
underlying counts and clustering of equations ([T]) andQ 
via the explicit forms of equations (11), ( jlsj ) and (201, 
and thereby produce the sample covariance, equation^ I, 
and Fisher matrix, equation (10). 



e.g. [57] and [55]). 

We subdivide the sky into 400 bins of 10 sq. degrees 
each, and calculate the counts and sample variance using 
mass bins of width log(AMobs) = 0.2 with the exception 
of the highest mass bin, which we extend to infinity. We 
set the width of our redshift bins to Az? = 0.1. These 
bin sizes imply 20 redshift bins and 5 bins of mass for 
the surveys 1 and 2. For Surveys 3 and 4, we divide the 
mass range lO^^.s < ^jopt ^ lO^^-^h'^MQ into 5 bins 
and use the same mass binning as the Surveys 1 and 2 
for M°P* > IO^^-^H-^Mq, with a total of 10 mass bins 
and 20 redshift bins. 

We assume fiducial cosmological parameters based on 
the fifth year data release of the Wilkinson Microwave 
Anisotropy Probe (WMAP5, [5S]). Thus, we set the 
baryon density, fli,h^ = 0.0227, the dark matter density, 
^mh^ = 0.1326, the normalization of the power spectrum 
at fc = 0.05Mpc"\ As = 4.625 x 10~^ the tih, n = 0.963, 
the optical depth to reionization, t — 0.087, the dark en- 
ergy density, J^de — 0.742, and the dark energy equation 
of state, w — —1. In this cosmology, erg — 0.796. 

With the exception of w, the cosmological parameters 
we used have been determined to an accuracy of a few 
percent. Extrapolating into the future, we assume 1% 
priors on all cosmological parameters except JIde and 
w. We used CMBfast [BD], version 4.5.1, to calculate the 
transfer functions. We do not explore time evolution of 
w in this work. In cosmologies with a time varying equa- 
tion of state, w(a), the relationship between halo model 
parameters and the error in principal components of u'(a) 
would differ from that reported here, with redshift terms 
having a larger impact. An important aim of next gen- 
eration surveys is to test w ^ —\, the degree to which 
dark energy differs from vacuum energy. 



D. Reference Model Surveys 



III. RESULTS 



We apply our tests to four distinct surveys consisting 
of a fiducial and three options. All are assumed to cover 
a sky area of 4000 square degrees and extend to a limit- 
ing redshift z^^ax = 2.0. The test surveys differ only in 
two parameters: the mass threshold, Mth, and the zero- 
redshift variance in the mass-observable relation, ctq. 

Our chosen survey parameters, given in Table [iTl rep- 
resent capabilities likely to be realized in the near future 
using sub-mm and optical/NIR observations. For exam- 
ple, the South Pole Telescope (SPT) is expected to de- 
tect clusters above Mth = lO^'*'^/i~^Af0 up to a redshift 
of 2 (see e.g. [5B]), with photometric redshifts available 
from DES+VISTA^ imaging, the Blanco Cosmology Sur- 
vey, and the Magellan Telescope. DES+VISTA wiU have 
internal capability to detect clusters optically with tech- 
niques known to work above Mth = lO^^'^hr^MQ (see 



http: //www. vista. ac .uk/ 



The baseline, absolute accuracy in f^DE and w mea- 
surements from the four test surveys are given in Table 
[TTJ The sharp priors columns represent the ideal condi- 
tion of perfect knowledge of the mass function and mass- 
observable relation, i.e., delta function priors on all the 
mass function and mass-observable relation parameters. 
The no priors columns give results assuming a high de- 
gree of ignorance in parameter values. 

From Table|lll we see that a factor 10*^'^ decrease in the 
mass threshold, Mth, improves the constraints on both 
JIde and w by a factor of ~ 4. The improvement results 
from the increase in cluster counts as well as an increase 
in the exposed range of the mass function. Surveys 2 and 
3 find a factor of 27 and 18 more clusters than the fidu- 
cial survey and survey 1, respectively. Increased scatter 
in Mobs results in an increase in counts because of the 
steepness of the mass function near the Mth. The nega- 
tive mass function slope implies that more objects will be 
up-scattered from below Afth than down-scattered. If the 
scatter is perfectly known, the increase in counts yields 



5 



TABLE II: Surveys parameters and constraints on cosmological parameters 

Sharp priors No priors 
Survey Mth[h^^MQ] ao Ntot cr(nDE) (^{w) o-(nDE) cr(w) 
Fid. lO^"'^ 0.2 8,400 0.010 0.050 0.91 2.19 

1 lO^''-^ 0.5 16,400 0.0083 0.039 0.82 1.81 

2 10^^-^ 0.2 359,600 0.0025 0.011 0.098 0.23 

3 10^^'^ 0.5 482,400 0.0023 0.0097 0.22 0.35 



better cosmological constraints [iTl [19] . This optimistic 
result must be interpreted with caution. A large mass- 
observable scatter may reflect poor selection, or contam- 
ination by projection or intrinsic sources. Projection is 
known to produce non-Gaussian features in the mass- 
observable relation [ST] that can bias constraints if not 
correctly accounted for (Erickson et al., in prep. |62)). 

The last columns of Table |ll] demonstrate that large 
degradations in cosmological accuracy result from essen- 
tially complete ignorance of the halo modeling parame- 
ters. We turn now to study the transition between the 
regimes of complete knowledge and complete ignorance 
by varying the prior uncertainties on the sixteen nuisance 
parameters. 



A. Degradation from Model Systematic Error 

For the parameters controlling the mass-observable re- 
lation (Afobs) and the mass-function/bias (MF/B), we in- 
troduce prior uncertainties, Cprior, that represent errors 
from previous observation or simulation. Because the 
Mobs and MF/B parameters have different dimensions, 
we define the priors on them differently, so as to make 
the prior uncertainties more directly comparable. 

We define the prior uncertainty on Mobs nuisance pa- 
rameters, cr^lioj: , such that the prior F^l-^^^ added to the 
i^^ diagonal element of the Fisher matrix is 



prior 



1 

prior 



(21) 



The prior uncertainty on an MF/B nuisance parameter, 
Xi, is defined as 



prior 



1 



^2 MF/B 
"prior , 



(22) 



With this definition, cr^^(^ corresponds to the prior frac- 
tional uncertainty on each mass-function/bias nuisance 
parameter. The uncertainties in the Mobs nuisance pa- 
rameters are already fractional since the mass-observable 
relation is defined in terms of the logarithm of the mass. 



The left plot of Fig. [T] gives a sense of the magnitude of 
the shifts in number counts as each of the Tinker MF pa- 
rameters are increased by a fractional amount, cr^^(.^ = 
0.1. We evaluate the mass function shifts for the fiducial 
cosmology at redshift z — 0.75, roughly the median red- 
shift of our model surveys. Solid lines show the effects of 
varying the constant terms while dashed lines vary the 
redshift-dependent factors. At lO^^ft.~^M0, the shift in 
number ranges from -1-0.2 (varying ag) to —0.3 (varying 
c). The derivatives are positive for the constant terms, 
with the exception of c, while the derivatives with respect 
to the z-dependent terms are negative. From the figure, 
we see that a 10% change in a causes very little change 
in the mass function. However, as we shall see in Sec. 
Ill B[ a and the other redshift evolution parameters {ux 



and Aj.) are nearly perfectly correlated (anti-correlated) 
with Que (w), suggesting that the redshift evolution of 
the mass function needs to be well understood to avoid 
degradations in cosmological parameter constraints. In 
the right plot, we show the dependence of the observed 
mass function on the Mobs nuisance parameters, at the 
same redshift as above. The impact of the Mobs nui- 
sance parameters increase with mass because the slope 
of the mass function increases with mass. The steeper 
the slope, the more significant is the imbalance between 
clusters up-scattered or down-scattered into a given mass 
bin due to bias or scatter in Afobs- 

Figs. [2] and |3] show contours of the multiplicative in- 
crease in the errors cr(i7DE) and <7{w), respectively, rel- 
ative to the baseline "Sharp prior" constraints given in 
Table |ll] The top left panel shows the fiducial survey, 
with survey 1 (upper right), 2 (lower left) and 3 (lower 
right) also shown. In all panels, contours show degrada- 
tion of the error in fluE or w by factors of 2-'/^, with j 
running from 1 to 8. 

The contours in Figs. [2] and [3] display similar shapes. 
Contours tend to intersect the axes at right angles be- 
cause of the very strong prior (10^'^) being imposed on 
one sub-space of parameters. The contour spacing gives 
the degradation in constraints as a function of prior on 
the complementary sub-space. For small systematic er- 
rors, the constraints on fluE and, especially, w degrade 
faster in the Mobs direction than in the MF/B direction. 
Since the Mobs systematic degrees of freedom only have 
redshift, not mass, dependence, this indicates a larger 
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1013 1014 IQ15 IQIS IQ14 IQ15 

halo mass [h"i Mg] halo mass [h^i Mg] 



FIG. 1: The sensitivity of the mass function to variation in Tinker mass function (left panel) and mass-observable (right) 
parameters is illustrated for redshift z — 0.75, roughly the median redshift of the surveys considered here. The change in 
the natural logarithm of number counts as a function of mass is shown while each of the labeled parameters is changed by a 
fractional amount, o-^f|^{^ = 0.1 (left panel) or a fixed amount, o-^l°^^? = ±0.1 (right panel). In the latter case, the bias terms 
are varied by 0.1 and the variance terms by —0.1. 



sensitivity to redshift evolution, particularly in the case 
of w. For example, consider a •\/2 increase in a{w) for 
the default survey (upper-left panel of Fig. [s]). For the 
case of perfect knowledge of the mass-observable priors, 
this degradation is reached when the mass function un- 
certainties are at the fractional level of ^ 50%. Contrast 
this to the case of perfect knowledge of the mass function 
priors, for which a ^/2 increase in aluu) occurs with only 
0.5% error in the mass-observable parameters. 

For the case of J7de shown in Fig. [2] the sensitivity to 
priors in Mobs and MF /B is more balanced for small val- 
ues of (Jpli^^ and (T^^{^ . For the fiducial survey, errors of 
^ 0.1 in the combined model parameters produce a factor 
2 increase in a{ilBE)- For a^l^^j: > 1.0, the sensitivity of 
the degradation to priors on Afobs parameters decreases 
sharply. This plateau reflects the ability of mass function 
shape and clustering information alone to jointly cali- 
brate Afobs parameters and JIde constraints. The same 
effect is noticeable in Fig. |3j though for w, the shape 
information provides the primary source of constraints 
of the Afobs parameters. When cr^^{^ is below the level 
of a few percent, the survey self-calibration is more ef- 
fective at constraining the A^obs nuisance parameters, so 
that increasing cr^li^^ above ~ 0.1 does not result in fur- 
ther degradation of the cosmological constraints. The 
corresponding plateau in the MF/B parameter direction 
is much less pronounced. Consequently, if no prior in- 
formation is available, mass-function/bias uncertainties 
dominate the error budget in cosmological parameters. 

The contours of fixed degradation shift as one considers 
the other surveys shown in the remaining panels of Figs. 



[2] and [3] Comparing the right panels (showing surveys 
with do — 0.5), to the left (surveys with ctq — 0.2), we see 
that increasing the default A/obs scatter tends to shift the 
contours inwards, signifying an increase in sensitivity to 
the priors. This increase in sensitivity offsets the smaller 
baseline error in JIde and w, leading to roughly constant 
errors in these parameters for fixed nuisance priors. The 
shift inwards is most noticeable for large values of Uprior 
for the surveys with Mth — 10^^-^h~^MQ. In the limit of 
flat priors on both MF/B and Afobs parameters, we see 
from Table n] that the increase in scatter is beneficial for 
the surveys with A^th = 10^'^'^h~^ Mq but detrimental 
to the surveys with A/th = 10^^'^h~^MQ. The former 
are dominated by shot noise, hence benefit from the in- 
crease in counts, but the latter are dominated by sample 
variance. 

Comparing the bottom panels (surveys with A/th — 
lQ^3-5fi~^MQ) with those above (surveys with A/th = 
lO^'*'^/i~^A/0), we see that the overall effect of decreasing 
the mass threshold is to increase the sensitivity to MF/B 
and Afobs priors. Both I^de and w degradation contours 
shift inwards by as much as an order of magnitude. The 
effect is most pronounced in the MF/B direction. For the 
case of w constraints with sharp Afobs priors, an interme- 
diate plateau emerges in the range cr^^{^ ~ 0.1 — 1. With 
lower Afth, there is more information in cluster surveys, 
which requires very accurate priors to fully extract. 

In summary, surveys with lower (better) baseline dark 
energy constraints have tighter requirements for priors 
on model systematic effects. For small cr^^or ^^'^ '''prkir^ 
(< 1), the larger degradation for tighter baselines nearly 
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FIG. 2: Plots of fractional degradation in a{QuE) for the Fiducial Survey {top left), Survey 1 {top right), Survey 2 {bottom left), 
and Survey 3 {bottom right). The contours correspond to regions where constraints on Qde are degraded by factors of \/2, 2, 
2^/2, 4, 4^/2, 8, 8\/2, 16, relative to the case of perfectly known nuisance parameters. See Table |ll] for the baseline constraints 

on f^DE- 



offset, in the sense that the absolute uncertainties in I^de 
and w at fixed c^l^^^ and cr^^{.^ are almost constant 
among the surveys. For high fiprior the better baseline 
does not compensate the more stringent prior require- 
ment because the surveys with larger scatter are more 
sensitive to the prior uncertainties. 

These results assume no mass-dependent evolution of 
the mass-observable relation. If we adopt the parameter- 
ization of [19' which consists of adding one parameter to 
describe the mass-evolution of the bias and three param- 
eters to characterize the mass-evolution of the scatter, 
we do not find significant qualitative changes. In sur- 
vey 3, for example, assuming cr^j^{.^ = ^-prior =0.1 the 



four additional parameters yield essentially no shift in 
the contours. In other regions of the Figure, results can 
be more noticeable, with shifts in the cr^^{^ or cr^l?^^ 
directions of up to a factor of 2. The rough shape of the 
curves is very similar. Degradations of cr{w) are some- 
what more sensitive to the inclusion of mass-dependent 
nuisance parameters. Assuming cr^^{^ = fpHor =0.1 
in Survey 3, addition of mass-dependence degrades con- 
straints by roughly \/2. The effects on the other surveys 
are even weaker since either they probe a smaller range 
in mass or they have smaller fiducial scatter. 
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FIG. 3: Plots of fractional degradation in a{w) for the Fiducial Survey {top left), Survey 1 {top right), Survey 2 {bottom left), 
and Survey 3 {bottom right). The contours correspond to regions where constraints on w are degraded by factors of \/2, 2, 2\/2, 
4, 4\/2, 8, 8\/2, 16, relative to the case of perfectly known nuisance parameters. See Table for the baseline constraints on w. 



B. Parameter correlations 

The halo modeling nuisance parameters have complex 
correlations among themselves and with the cosmological 
parameters. While one could imagine adopting a more 
orthogonal parameterization with potentially fewer pa- 
rameters, the interpretation of the results would be hard 
to relate to the functional forms currently in use. Thus, 
we pursue the direct approach of exploring the Fisher 
correlations of the full parameter set. We begin by iso- 
lating the Mobs and the MF/B sub-spaces separately, us- 
ing sharp priors on the complementary sub-space, then 
consider the full parameter covariance matrix under the 
assumption of a 10% prior uncertainty on the halo mod- 



eling parameters. 

First, we add sharp priors to the MF/B nuisance pa- 
rameters and consider the correlations between the Mobs 
nuisance parameters and the dark energy parameters, 
shown in Fig. [4] (top left). The labels of the Mobs and cos- 
mological parameters occupy the diagonal of the correla- 
tion matrix. As described in Sec. |II D| the cosmological 
parameters other than f^DE and w have 1% priors, and 
this level is small enough to remove almost all of their 
correlations with the other parameters. The exception 
is the normalization of the primordial power spectrum, 
\n{As), which strongly correlates with ilpE despite the 
1% prior. 

From the plot, we see that the constant nuisance pa- 
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FIG. 4: Correlations between cosmological parameters and {top left) M^bs nuisance parameters, and {top right) MF/B nuisance 
parameters. The nuisance parameters not shown in each plot are fixed by infinitely sharp priors. Cosmological parameters, with 
the exception of S^de and w have 1% priors and are not shown unless significant correlation with other parameters is present. 
The bottom plot shows the correlation between Mobs, MF/B, ^de, and w. We have added priors to all nuisance parameters 



such that a„ 



MF/B 



0.1, which very roughly corresponds to our present knowledge of these parameters. 
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rameter of the mass bias (Bq) is strongly correlated 
with r^DE whereas the redshift-related parameter (Bi) 
is markedly anti-correlated with J7de- The variance nui- 
sance parameters (ctq, 5*1, 5*2 , S3) do not correlate as 
strongly with fluE- In contrast, the equation of state w 
is mostly correlated with the redshift-dependent variance 
parameters, and somewhat correlated to Bq. These cor- 
relations are consistent with the trends seen in Fig. 15 of 
[19j . showing the dependence of riDE and w constraints 
on the priors on different Mobs parameters. 

On the {top right) plot of Fig. |4]we examine the cor- 
relations between the mass function nuisance parameters 
and the cosmological parameters. As in the above case, 
we apply infinitely sharp priors to all nuisance parame- 
ters not shown, and only plot the parameters with signif- 
icant correlations. In this scenario, the ln{As) correlates 
more weakly with all other parameters. The redshift- 
dependent parameters (a, a^, A^) show strong positive 
correlations with fluE, and strong negative correlations 
with w. The exponential cutoff parameter c is also notice- 
ably correlated with w and ^Ibe- The redshift-dependent 
nuisance parameters are also strongly correlated among 
themselves, as are the constant parameters (Aq, ao, and 
60)) which are very weakly correlated with all other pa- 
rameters. 

In the bottom panel of Fig. |4j we investigate the full 
correlation matrix after imposing priors corresponding to 



Moi 



MF/B 



VroT = <^prior = 0.1 to the Mobs and MF/B parame- 
ters (along with the 1% cosmological priors). The ln(v4s) 
correlations are very weak, hence, we do not show them. 
The choice of priors on the halo modeling parameters 
is admittedly crude, based roughly on the present-day 
understanding of the mass-function, bias, and the mass- 
observable relation. We caution that the correlations are 
a strong function of the imposed priors, so the chosen 
case is illustrative rather than definitive. 

With these priors, the uncertainties in the Mobs nui- 
sance parameters dominate the error budget in I^de and 
w. This result is consistent with the shape of the contours 
near (0.1, 0.1) in Figs. 2]and[3] Thus, the correlations 
between the Mobs and dark energy parameters are more 
pronounced, and, in general, resemble the correlations, 
performed under sharp MF/B priors, displayed in the 
top left panel. An exception is the correlation between 
the mass bias constant, Bq, and I^de, for which the cor- 
relation in the full treatment is substantially weaker. 

The mass function parameter correlations are substan- 
tially different from their isolated treatment. In partic- 
ular, the redshift evolution parameters {Ar^, a^, and a) 
largely disappear, and only the constant parameters (Aq, 
ao and 60) contribute appreciably to the dark energy er- 
ror budget. That this behavior differs from the isolated 
case is not too surprising. In the isolated case, the Mobs 
parameters are assumed to be perfectly known, so the 
only redshift evolution remaining in the model to com- 
pete with dark energy is contained in the MF/B param- 
eters. In the full case, the assumed prior level of Mobs 
uncertainty is sufficient to dominate the evolutionary be- 



havior, leaving the primary shape parameters of the MF 
as the means by which this sector affects dark energy 
constraints. 

The full analysis also includes halo bias nuisance pa- 
rameters. The parameters 5c and abias strongly correlate 
with each other and also display weakly negative correla- 
tions with both r^DE and w. The parameter pc is virtually 
uncorrelated with all the other parameters. 

As we have mentioned before, the detailed behaviors 
are a consequence of the priors applied to the nuisance 
parameters. Had we used flat (very weak priors), the 
MF/B nuisance parameters would dominate the error 
budget, since self-calibration provides some constraints 
on the Mobs parameters. Then, the correlations of the 
mass function parameters between themselves and the 
cosmological parameters would more closely resemble 
those seen at the top right panel. The bias parameters 
would also exhibit somewhat stronger correlations. 

To further illustrate this sensitivity to priors on nui- 
sance parameters, we explore the contributions of dif- 
ferent sets of MF/B and Mobs nuisance parameters to 
degradations in cr(riDE) and a{w). The top row of Fig. 
[5] shows the fractional degradation of a{^EiE) and a{w) 
with respect to the baseline of the fiducial survey for the 
cases where (1) all the MF/B parameters are allowed 
to vary, (2) only the MF parameters are free, and (3) 
only the redshift-independent MF parameters are free. 
As before, wc fix the priors on Afobs parameters so that 

Mobs ^0.1. 
prior v^.-L. 

In the upper left plot of Fig. [5] we see that the redshift 
evolution MF nuisance parameters {A^^ a^, and a) dom- 
inate the degradation of a{^luE) for moderate to large 
error, 0.1 < cr^^{.^ < 1- Below cr^^(^ ~ 0.1, the con- 
stant MF parameters (^Oi bo and c) are the most 
relevant. The bias parameters are relevant in the high- 
uncertainty range, 1 < cr^JJ^^ < 10. When cr^^(^ = 0.1, 
the parameters not related to redshift evolution, i.e., Aq, 
ao, and bo dominate, as seen in Fig. [4] In the upper 
right plot, the constant MF nuisance parameters dom- 
inate the w constraints up to cr^^^^ ^ 1. For larger 



MF/B 
prior ' 



the redshift evolution parameters become more 
important. The bias parameters are moderately relevant 
in the range cr^i^{^ ^ 0.03 — 1. 

The bottom row of Fig. [5] shows the degradation of 
51de and w constraints as we vary priors on Mobs param- 
eters while keeping cr^^^^ fixed at 0.1. From the lower 
left plot, we see that the parameters related to the red- 
shift evolution of the mass-observable relation are the 
most important for virtually the entire interval we exam- 
ine. In particular, the redshift evolution of the mass bias 
is the most relevant for JIde constraints. In the lower 
right plot we see that the redshift evolution of the mass 
variance dominates w constraints for small uncertainties 
in the priors. For a^l^^j: > 0.1, the constant part of the 
variance, ctq, as well as the redshift evolution of the bias, 
dominate. Constraints are almost independent of priors 
on the constant bias term. Bo, but this term does affect 
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the power spectrum normalization, hi(As). 



IV. DISCUSSION 

The results presented in this paper make a variety of 
assumptions of various degrees of relevance, which must 
be interpreted with caution. In this section we discuss 
the generality of some of our assumptions. 

Parameterization of the mass-observable relation We 
parameterized the redshift evolution of the variance in 
Mobs using a cubic polynomial. Lima & Hu (2005) [T7] 
show that the cubic polynomial is almost as complete a 
description as having fully independent scatter in ^ 20 
redshift bins. Hence, we feel that our parameterization 
is conservative with regards to redshift evolution. There 
is no physical motivation for this choice, however, and 
if simpler parameterizations describe the data well, then 
constraints would improve - and the sensitivity to the 
uncertainty in the Mobs parameters would decrease. As 
a test, we eliminated the quadratic and cubic terms in 
the scatter in Mobs- The contour lines of Figs. [2] and [3] 
shifted upwards by as much as factors of 5 in cTpj-ior- 

We did not include mass evolution of the mass- 
observable relation since evidence supporting this as- 
sumption from observations and simulations is currently 
weak. For surveys with high Afthj we have checked 
that cosmological constraints are virtually unaffected by 
adding a cubic evolution of the mass scatter plus a linear 
evolution of the mass bias. For surveys with low Mth, the 
mass terms cause an increase in sensitivity to both A/obs 
and MF/B parameters, which results in a shifting of the 
contours of Figs. [2] and [3] inwards. The additional Mobs 
nuisance parameters correlate strongly with the MF/B 
nuisance parameters, so that the sensitivity of both sets 
of parameters to the priors vary in similar fashion. 

Our results are based on a quite generic procedure for 
adding priors. We did not strive for optimal dark en- 
ergy constraints, or to accurately reproduce what spe- 
cific observations and simulations might return. Under- 
standing the optimal priors needed - for fixed observa- 
tional/simulational costs - can be very valuable |63j . and 
such detailed studies are important complements to the 
generic treatment presented here. 

As mentioned previously, we did not consider varia- 
tions in w in this paper. The sensitivity of the Dark En- 
ergy constraints to the redshift-dependent nuisance pa- 
rameters should increase, though it is not obvious that 
the relative importance of the Mobs and MF/B nuisance 
parameters would change significantly. The subject of 
cluster constraints on time-varying w, parameterized as 
Wfj/wa, is addressed in papers such as Wu et al. |51j and 
Cunha et al . 



V. CONCLUSIONS 

We investigate the sensitivity of dark energy con- 
straints from clusters of galaxies to halo modeling un- 
certainties in the mass function, clustering bias, and 
the mass-observable relation. We find that mass- 
observable uncertainties dominate the error budget for 
both f2DE and w constraints for surveys with higher 
mass-thresholds, such as SZ surveys, assuming prior un- 
certainties of order 0.1 on both mass-observable and 
mass-function/bias (MF/B) nuisance parameters. For 
surveys with lower mass-thresholds, the uncertainties 
in the mass-observable and MF/B nuisance parameters 
are more comparable, depending on the degree of prior 
knowledge assumed. 

The variations in the sensitivity to the prior uncer- 
tainties are offset by the different baseline constraints of 
each survey. Not surprisingly, surveys with lower (better) 
baseline constraints have tighter requirements for priors 
on model systematic effects. 

We examine the correlations between the nuisance and 
cosmological parameters for the fiducial survey for sev- 
eral different assumptions about prior knowledge of the 
nuisance parameters. If the mass-observable relation 
is perfectly known, the mass function parameters show 
strong positive (negative) correlations with JIde and w. 
If the mass function parameters are known to ~ 0.1, then 
"■(J^de) and a{w) are dominated by the constant param- 
eters of the mass function (^qi ^Oi and 69). When the 
mass-function is held fixed, JIde is most sensitive to the 
normalization (bias) of the mass-observable relation and 
its redshift evolution, whereas w is more sensitive to the 
redshift evolution of the variance. 

We only consider individually self-calibrated cluster 
surveys. But as [11119) show, cross-calibration is a power- 
ful tool to improve knowledge of the mass-observable nui- 
sance parameters and thereby tighten cosmological con- 
straints. Strategies to optimize follow-up observations of 
cluster samples Wu et al. are relevant in this regard. 
The effects of using cross-calibrated cluster surveys would 
be to decrease the sensitivity to the mass-observable pa- 
rameters, thereby increasing the relative sensitivity to 
the MF/B parameters. Our results adopted a specific 
parameterization of the mass function based on the as- 
sumption of collisionless dynamics of dark matter, with 
no gas physics. Baryonic effects can be large on the scale 
of clusters Stanek et al. j47l , and simulations exploring a 
broader range of baryon physics behavior are needed to 
ensure that the uncertainties in the mass function do not 
play a more significant role in limiting our knowledge of 
dark energy. 

In conclusion, our results further illuminate challenges 
to precision cosmology with galaxy clusters. Observers 
and simulators should focus on characterizing the form 
of the mass-observable relation, since having the correct 
parameterization is essential to avoid biases in derived 
cosmological constraints. While we find that a mass func- 
tion accuracy of about 10% is a sub-dominant source 
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FIG. 5: Degradation of constraints on (left) Qde and (right) equation of state ui as a functi on of the prior on the uncertainty 

Mo 



in MF/B and Mobs nuisance parameters: a. 



MP/E 



and 



respectively (defined in Eqs. 



21 



and 



22 1. Plots in the top row 



assume prior uncertainty on the A/obs parameters of o-p^.P^j.'' — 0.1. The solid black lines assume no priors on any of the MF/B 
nuisance parameters. The dashed red lines assume sharp priors on the three bias parameters, and the solid green lines assume 
sharp priors on the redshift evolution mass function nuisance parameters (Ax, a^, and a). In the range of (T^^{^ plotted, the 
green line is unaffected if priors are applied to the bias parameters or not. Plots in the bottom row assume prior uncertainty 
on the MF/B parameters of o-^l°^^ — 0.1. The solid black lines assume no priors on any of the Mobs nuisance parameters. The 
dashed red lines assume sharp priors on the three parameters describing the redshift evolution of the mass variance. The dotted 
blue lines assume sharp priors on all four mass variance nuisance parameters (ctq, Si, S2, and S3), and the solid green lines 
assume sharp priors on all Afobs nuisance parameters except the constant bias term Bq. 
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of degradation of cosmological constraints, it still con- 
tributes measurably, particularly for surveys where sig- 
nificant cross-calibration or targeted follow-up will be 
possible. 
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